Unified Chew-Mandelstam SAID analysis of pion photoproduction data 
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hadroproduction data, has been achieved in a Chew-Mandelstam parametrization which is consistent 
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I. INTRODUCTION 

A wealth of 7-/V — > ttN data, for single- and double- 
polarization observables, is anticipated from electromag- 
netic facilities worldwide over the coming months and 
years. These data will be pivotal in determining the 
underlying amplitudes in complete experiments, and in 
discerning between various microscopic models of multi- 
channel reaction theory. 

The focus of precision electromagnetic measurements, 
over the nucleon resonance region, is to more fully map 
the non-perturbative regime of quantum chromodynam- 
ics, the fundamental theory of the strong interaction, to 
shed light on its confining and chiral symmetry breaking 
properties. These electromagnetic data take the field to 
the next and necessary level of precision. This is required 
in order to obtain a theoretical description of the nucleon 
that both explains and subsumes the simple constituent 
quark model, which has provided a qualitative picture 
of nucleon structure and reactions. The expected data 
heralds an era of precision hadron spectroscopy, partic- 
ularly for baryons, and has ushered in a renaissance in 
hadronic reaction theory. Significant refinements in the 
quality and quantity of available data offer the oppor- 
tunity to develop more sophisticated models of hadronic 
reactions, constrained by fundamental principles of field 
theory, such as unitarity and gauge invariance, which 
have model dependencies under better control, if not 
eliminated. Such a complete and successful phenomenol- 
ogy would appear to be a prerequisite for a deeper un- 
derstanding in terms of quarks and gluons [l[ . 

The present manuscript details multipole analy- 
ses of the single-pion photoproduction data using a 
parametrization form related to, but an improvement 
upon, previous SAID parametrizations [2tJ]- The 
energy-dependent (ED) analysis is performed over the 
center-of-mass energy (W) range from the near-threshold 
region to about 2.5 GeV, including resonances through 



the fourth resonance region. We also generate single en- 
ergy solutions (SES), which fit the data over narrow en- 
ergy bins assuming phase information obtained from the 
ED solutions. The relations between ED and SES fits 
have been extensively studied in Ref. |5| . A detailed dis- 
cussion of amplitude and observable conventions is also 
given in this source. 

The fitted pion-photoproduction database is identical 
to that used in our most recent [|| SN11 analysis, based 
on the standard SAID parametrization. In the following 
section, we compare the previous and present SAID fit 
forms used to analyze these data. Extracted multipoles 
are compared to previous SAID fits, and those from other 
groups, in Sec. IIII1 Our results and their implications are 
summarized in Sec. IIVI 



II. FORMALISM 

The Chew-Mandelstam (CM) energy-dependent (ED) 
parametrization for the hadronic T matrix, described in 
Ref. 0], has been used in a recent coupled-channel fit of 
irN elastic scattering and ttA^ — s> rjN reaction data. It 
gives a realistic description of the data with x 2 per datum 
better than any other parametrization or model, to our 
knowledge 0, H[ . The parametrization form used in this 
fit is given as 



(1) 
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where a, /3, and a are indices for the considered chan- 
nels, 7riV, 7rA, pN , and wN . This parametrization has 
been discussed in Refs. 0, Hcj |. Given the success of 
this approach in the hadronic two-body sector, its ap- 
plication to the study of meson photoproduction is war- 
ranted. The main result of the present study is the use 
of the information encoded in Eq. (p} by employing the 
factor [1 — K(W)C(W)}~ 1 (called the "hadronic rescat- 
tering matrix") in the photoproduction parametrization 
form. 

The CM form of Eq. |T]) may be extended to include 
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the electromagnetic channel as: 

T ai = Y^-KC\^K ai . (2) 

a 

Here, 7 denotes the electromagnetic channel, jN, and 
a denotes the hadronic channels which appear in the 
parametrization of the hadronic rescattering matrix, 
[1 — KC}^ 1 . Note that by sharing the common factor, 
[1 — KC]^ 1 which encodes, at least qualitatively speak- 
ing, the hadronic channel coupling (or rescattering) ef- 
fects, Eqs. (HJ) and constitute a unified approach to 
the problem of parametrizing the hadronic scattering and 
photoproduction amplitudes. 

We pause here to make several remarks about the an- 
alytic form of the parametrization and its use in the 
present study. We first note that since the CM K-matrix 
K ai {yV) is a polynomial in the center-of-mass energy, W 7 
an entire function, non-analytic points in the complex- W 
plane are all a result of the hadronic rescattering ma- 
trix, [1 — KC}^ 1 . This matrix has branch points and 
poles consistent with two-body and quasi-two-body uni- 
tarity 1 1 1 . The quasi-two-body channels, nA and pN 
model the three-body ttttN channel only approximately. 
We use Eqs. ((TJ and © as follows. The parameters of 
the hadronic CM K-matrix, K aj p, where a and /3 may in- 
clude (depending on the partial wave) nN,TrA,pN, and 
r)N are fixed by fitting the irN — > irN and ttN —> r]N 
data as in Ref. [6j. The 7r-photoproduction data is then 
fitted [IH by varying only the parameters of the electro- 
magnetic CM K-matrix elements, K aj , where a includes 
the channels ttN, ttA, pN, and r/N. 

This approach differs markedly from that adopted in 
Refs. @, II [H, d. There, the fit form, motivated by a 
multichannel Heitler K-matrix approach [15| . 

M = (Bom + A)(l + iT 7rN ) + BT 7rN , (3) 

was modified to include a term 

{C + iD){lmT^ N -\T^ N \ 2 ), (4) 

where T^n is the elastic TtN scattering partial-wave am- 
plitude associated with the pion-photoproduction multi- 
pole amplitude M . The added piece, which grows with 
the ttN reaction cross section, was found to improve the 
fit at energies above the rjN threshold. Each of the phe- 
nomenological terms A to D was parameterized as a poly- 
nomial in energy, having the correct threshold behavior. 

In the new form, terms A and B have been effectively 
replaced by a single CM K-matrix element, K %y . Con- 
versely, the influence of channels opening above ttN is 
now (more correctly) associated with individual channels 
(i]N, 7rA, and pN) rather than a single term. Large can- 
cellations found to occur between the Born, A, and B 
terms suggested that a more economical parametrization 
is possible 15]. In fact, the CM K-matrix form provides 
a better overall fit to the data with fewer free parameters, 
as we will show in the next section. 



III. FIT RESULTS AND MULTIPOLE 
AMPLITUDES 

In Table HI the fit quality and number of searched 
parameters is compared for the two fit forms discussed 
in the previous section. Here we have used the same 
database to 2.7 GeV as was used in generating solution 
SN11 Q. The present CM12 form requires fewer param- 
eters to achieve a slightly better data fit. The energy- 
dependence of this result was tested by repeating the 
analyses over three different energy ranges. 

As the fit employs polynomial functions for A to D, of 
SN11, or the electromagnetic CM K-matrix elements of 
CM12, a subjective criteria is required to determine the 
order of polynomials fitted. In the fits to 2.7 GeV, the 
order of polynomial functions was increased until further 
additions improved the overall x 2 by 50 or less. This 
same criteria was used in both the SN11 and CM12 fits, 
in order to have a basis for comparison. As more pa- 
rameters were searched, their ability to improve the fit 
diminished. In the fits to lower energy limits, parame- 
ters were removed in steps, again with the condition that 
removing a parameter should not increase the overall x 2 
by more than 50. 

In Figures HHS we compare SN11 and CM 12 to 
MAID07 [HI and Bonn-Gatchina [13 fits. Only proton- 
target multipoles are presented, as changes in the 
neutron-target database are likely to alter these fits in 
the near future. For resonances with a canonical Breit- 
Wigner shape, such as the A(1232)P 33 , A(1520)L>i 3 , 
7V(1680)Fi 5 , and A(1950)F 37 , all solutions agree fairly 

well in the neighborhood of the resonance energy. The 

1/2 

fit CM12 deviates significantly from SN11 in the E ' + 
multipole. The CM12 phase behavior, from threshold up 
to the peak of the N(1535)Sn resonance, differs from 
SN11 and MAID07, but is qualitatively similar to the 
Bonn-Gatchina result, as shown in Figure [7] While 
both SN11 and MAID07 essentially follow the Su pion- 
nucleon phase up to the r]N threshold cusp, the CM12 
and Bonn-Gatchina fits depart from this phase above the 
two-pion production threshold. 

Some structures occuring in SN11, between threshold 

and the first resonance energies, are missing or dimin- 

3/2 

ishcd in CM12. Examples are the real parts of M x _ , 

1 /2 1/2 

Mjl , and E 3 _ . In each of these cases, the CM12 fit 
more closely resembles the MAID07 and Bonn-Gatchina 
results. This reflects the fact that the replacement of the 
phenomenological terms A and B of Eq. ([3]) with a single 
CM K-matrix element in Eq. ^ is more form-restrictive. 

In Fig. we compare the % 2 /data of the ED and SES 
fits over energy bins used to generate the SES. As in 
Ref. we see a noticeable increase in the x 2 differ- 
ence above about 800 MeV. In Fig. [3J we compare x 2 
values for SES generated from the SN11 and CM12 fits 
over identical energy bins. As multipole phases are fixed 
for multipoles searched in generating these SES this 
serves as a comparison of the, often quite different, phases 
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found in SN11 and CM12. While the SN11 SES achieve 
a better fit in the near-threshold region, between 700 and 
1100 MeV the CM12 SES consistently give the better fit. 

TABLE I. xVdata and number of searched parameters (N p ) 
compared for fits to pion photoproduction data over varied 
energy ranges. The fit form used for solution SN11 [J is 
compared to the Chew-Mandelstam form, CM12. See text 
for details. 



Solution 


Energy limit 
(MeV) 


X 2 / N Data 


N p 


SN11 


2700 


2.08 


209 


CM12 


2700 


2.01 


200 


SNlla 


2100 


1.96 


206 


CM12a 


2100 


1.88 


194 


SNllb 


1200 


1.69 


175 


CM12b 


1200 


1.64 


166 



IV. SUMMARY AND CONCLUSION 

We have fitted the single-pion photoproduction 
database utilizing a parametrization consistent with the 
Chew-Mandelstam form used in our previous fits to irN 
scattering and r)N production data. This new fit has a 
number of interesting features. It is more economical, 
using fewer parameters to obtain a slightly better over- 
all fit. Some low-energy structures, not seen by other 
groups, have disappeared in the present fit. The phase 

1/2 

behavior of the E 0+ has changed significantly and now 
is qualitatively similar to the Bonn-Gatchina result. 

Comparison of the ED and SES fits shows, as was 
found in Ref. Q, a rise in % 2 difference, evaluated over 
narrow energy bins, above about 800 MeV in the pho- 
ton energy. This could be related to the energy limit of 
MAMI-B, which has contributed a significant fraction of 
the precise data below 800 MeV (thus, a data issue) or 
due to the treatment of channels above single-pion pro- 
duction (a model issue). The comparison of SES fits, 
derived from the ED SN11 and CM12 solutions, suggests 
that the CM12 multipole phases, held fixed in SES, are 
prefered in the intermediate-energy region. 

In Table HH we compare photo-decay couplings ex- 
tracted from CM12 and SN11, using the method of 
Ref. [4J, to the average of PDG values. As expected, 
the V(1535)»Sii shows a large increase to 128±4 (in 
GeV)~ 1/2 x 10~ 3 units), compared to 105±10 from the 
Bonn-Gatchina group, and 118 found in Eta- MAID [20| . 
As in the SN11 analysis, the V(1650)Sii is very difficult 
to fit using this procedure. The fit prefers a larger Breit- 
Wigner mass and width, and a lower value for T^/T. Al- 
lowing for these variations, as reflected in its uncertainty, 
the value quoted here should be taken only as a rough 
estimate. 

Other couplings are generally either close to values 



TABLE II. Resonance parameters for N* and A* states from 
the SAID fit to the nN data f(| (second column) and pro- 
ton helicity amplitudes A 1/2 and A 3 / 2 (in [(GeV)~ 1/2 x 10~ 3 ] 
units) from the CM12 solution (first row), the SN11 [3| solu- 
tion (second row), and average values from the PDG10 plj 
(third row). 



Resonance 


nN SAID 


A 1/2 


A3/2 


Af(1535)Sn 


W fl =1547 MeV 
T=188 MeV 
lVr=0.36 


128±4 
99±2 
90P30 




Ar(1650)Sii 


Wj?=1635 MeV 
T=115 MeV 

r^/r=i.oo 


55±30 
65±25 
53±16 




iV(1440)Pn 


W R =U85 MeV 
T=284 MeV 
lVr=0.79 


-56±1 
-58±1 
-65±4 




AT(1720)Pi 3 


W fl =1764 MeV 
T=210 MeV 

r,/r=o.o9 


95±2 
99±3 
18±30 


-48±2 
-43±2 
-19±20 


Ar(1520)Di3 


Wfl=1515 MeV 
T=104 MeV 

r^/r=o.63 


-19±2 
-16±2 
-24±9 


141±2 
156±2 
166±5 


Af(1675)Di 5 


W fl =1674 MeV 
T=147 MeV 
lVr=0.39 


13±1 
13±2 
19±8 


16±1 
19±2 
15±9 


iV(1680)Pi5 


W fl =1680 MeV 
T=128 MeV 
IVL=0.70 


-7±2 
-13±3 
-15±6 


140±2 
141±3 
133±12 


A(1620)S 3 i 


W fl =1615 MeV 
T=147 MeV 

r^/r=o.32 


29±3 
64±2 
27±11 




A(1232)P 33 


W R =1233 MeV 
r=119 MeV 

r^/r=i.oo 


-139±2 
-138±3 
-135±6 


-262±3 
-259±5 
-250±8 


A(1700)Z> 33 


W fl =1695 MeV 
T=376 MeV 

r, r /r=o.i6 


105±5 
109±4 
104±15 


92±4 
84±2 
85±22 


A(1905)P 35 


W fl =1858 MeV 
T=321 MeV 

iVr=o.i2 


19±2 
9±3 
26±11 


-38±4 
-46±3 
-45±20 


A(1950)P 37 


VTr=1921 MeV 
T=271 MeV 
T,r/T=0.47 


~83±4 
-71±2 
-76±12 


-96±4 
-92±2 
-97±10 



found in SN11, or are within the PDG ranges. The 
V(1520)P>i3 A 3/2 (141±2) has dropped below the PDG 
range (166±5), but is above the Bonn-Gatchina value 
of 131±10. The V(1720)Pi 3 couplings remain very un- 
certain, as is clear from the multipole plot comparison. 
Here the present PDG range is certainly too narrow - 
the Bonn-Gatchina result for is 150±30, compared 
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FIG. 1. (Color online) 1=3/2 multipole amplitudes from threshold to W = 2.43 GeV {E 1 = 2.7 GeV). Solid (dash-dotted) 
lines correspond to the SN11 (CM12) solution. Short-dashed (dashed) lines give BG2010-02 solution P3 (MAID07 [l]|, which 
terminates at W=2 GeV). Vertical arrows indicate resonance energies, Wr, and horizontal bars show full (F) and partial (T^n) 
widths associated with the SAID -k N solution SP06 [(| . 
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to the PDG range -19 ± 20. For A 1/2 , the PDG range 
is 18 ± 30, with the CM12, SN11, and Bonn-Gatchina 
results all near 100. 
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FIG. 2. (Color online) Notation of the multipoles is the same as in Fig. [T] 
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FIG. 3. (Color online) Notation of the multipoles is the same as in Fig. [T] 
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FIG. 4. (Color online) Proton multipole 1=1/2 amplitudes from threshold to W = 2.43 GeV (E-y = 2.7 GeV). Notation of 
the solutions is the same as in Fig. [T] Vertical arrows indicate resonance energies, Wr, and horizontal bars show full (r) and 
partial (r^-jv) widths associated with the SAID irN solution SP06 [a]. 
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FIG. 5. (Color online) Notation of the multipoles is the same as in Fig. [4] 
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FIG. 6. (Color online) Notation of the multipoles is the same as in Fig. [4] 
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of the solutions is the same as in Fig. [1] The thick solid line 
corresponds to the SAID ttN solution SP06 [6j for the Su 
phase shift. 
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